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Introduction 

Engineering systems are often assembled from a number of independent 
components which must function harmoniously according to the desired 
performance of the overall system. Although each component may function 
satisfactorily by itself, an assembly of such components may fail to meet 
the desired overall requirements. It is therefore necessary to examine 
such an assembly of components from the point of view of the system 
objectives. 

System engineering is an approach which recognizes the many variables 
of the system and interrelates them to the overall function of the system. 
The objectives of the system may be diverse aims such as low cost, reliability, 
vibration- free operation, or any other requirement to be specified. To attain 
the system objectives, modification of the subf unctions of the components may 
be required. If the system is to be optimized for a given performance, the 
effect of variation of each of the variables in the system must be examined. 
Of frequent concern is the question of system stability. Thus, fundamental 
aspects of many fields may be involved in a systems study* 
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The Pogo Problem 

One of the interesting examples of systems instability was that 
encountered by the behavior of the Titan-2 missile shortly after 
lift-off. The Titan-2 missile is a liquid propellant rocket 103 ft. 
tall and weighing 165 tons with a first stage thrust of 430,000 lbs. 
It was designed as a military weapon to carry a nuclear warhead 6,000 
miles or more. 

When the Gemini program was conceived, the Manned Spacecraft 
Center decided to use the Titan-2 missile to boost the two-man teams 
of astronauts into orbit. However, serious vibration problems had to 
be overcome for the booster. After about 100 seconds of flight, severe 
longitudinal oscillations, called the M Pogo bounce 11 between 11 and 14 
c.p,s., with accelerations up to 5 g's, would develop in the Titan-2, 
For an unmanned flight such vibrations could be tolerated,, However, 
to an astronaut, such vibration levels speiled trouble > For example, 
it is known that vibration levels of only 1/2 g at frequencies near 
10 c«p.s. would seriously impair one's vision, so that an astronaut 
subjected to such environment could not read the instrument panel. 
Moreover, such a ride would be a most unpleasant physiological experience* 

This longitudinal oscillation was known to exist in all liquid 
propellant missiles but its cause and cure were unknown in 1962. Various 
fixes were proposed and tried, and in one test a 15 million dollar missile 
was destroyed in flight due to excessive vibrations. After an intensive 
study program costing 3 million dollars, the troublesome f, Pogo bounce 11 
was identified as coupled oscillation between the structure and the 
rocket-propellant-combustion system. Components when tested separately 
would behave according to expectation and the difficulties experienced 
in operation as a system became predictable only by the systems approach 
of analysis. As in most cases, understanding the cause leads to a cure, 
and the "Pogo bounce 11 in the Titan-2 booster was successfully reduced to 
levels acceptable to the Gemini program. 

Systems Analysis of the Liquid Propellant Missile 

Figure 1 illustrates the essential components of the Titan-2 
propulsion system with only one of the two engines shown. The mono- 
propellant system shown in Figure 2 is a simplification which will 
adequately describe the instability known as the Pogo-bounce and 
illustrates the systems approach to be used. The components of the 
system are; 

(a) 



(b) 
(c) 

(d) 

(e) 
(f) 

2 



Propellant tank with its liquid propellant 
(We will refer to the propellant as fuel, 
although it could stand for the oxidizer as well.) 
Fuel feedlline 
Turbo- pump 

Discharge line and injector 
Rocket engine 
Missile structure 
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To simplify the analysis we will lump the mass of the system at 
two points where the displacements Y- and Y are measured from an 
inertial reference, (assumed to be the Earth). Thus, will include 
the mass of the turbo-pump and the rocket engine, while M? will represent 
the variable mass of the fuel in the tank. The two masses are connected 
by a spring and damper representative of the intervening structure, and 
the system is to be adjusted to result in the same natural frequencies 
as the first and second modes of the missile. 

The behavior of the system can be briefly described in a qualitative 
manner as follows, A disturbance in the fuel feed system is transmitted 
to the injector of the rocket in the form of a pressure variation and a 
corresponding variation in the fuel flow rate to the combustion chamber , 
The resulting fluctuation in the thrust introduces a variation in the 
structural excitation which causes relative motion of the fuel tank with 
respect to the feed line, the turbo-pump and the engine. With proper 
phase and period of vibration, the sequence of events described may 
continue to increase in amplitude to proportions unacceptable to. the 
design , 

Pertirient to.the problem is the stability of the system performance 
and how it may be altered by design modifications. The systems approach 
to be used here is that of establishing the transfer function of each 
component and assembling these into a block diagram circuit representative 
of the overall system. Such a circuit may then be analyzed by the various 
techniques available in the treatment of closed-loop systems. 

Component Transfer Function 

It is evident that there are several variables to be considered in 
the system. They can be identified by isolating each component and 
examining its input and output relationship which will be expressed in 
terms of a transfer function. For its development we will consider the 
fluctuation in the quantity under consideration to be a small perturbation 
from its steady state value and linearize the analysis by ignoring second 
order terms of the perturbations * The transfer function is then defined 
by the Laplace transform of the ratio of the output to input with zero 
initial cond it ion „ 

(a) Fuel Feed System 

The flow rate in the fuel feed system is affected by the acceleration 
of the structure and the pressure variation in the fuel. Referring to Fig, 2, 
the tank and its fuel are assumed to have acceleration Y2 whereas the fuel 
feed line is assumed to have acceleration Y\ equal to that of the pump and 
engine. Sections of this system from the tank to the suction inlet of the 
pump will be first analyzed by using Newton's second law, Bernoulli ! s 
equation, and the continuity principle of mass flow rate. Capital letters 
will be used to represent the total quantity, which will be separated into 
its steady state value with an overhead bar and its fluctuation' indicated 
by the lower case letter (i,e c P » P + p) > 
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Between points (3) and (2) Newtons equation of motion is 

where 

= pressure —lbs/ft 2 

3 

a weight density of fuel-^lbs/ft 
« head -vf t 

2 

* acceleration of gravity "^32 . 2 ft/sec 
» absolute acceleration of tank and fqel-^ft/sec 




For the stea dy st ate condition the missile is accelerating with t 
acceleration Y .under steady pressures P: and'iwe sim{>ly.:ob£ain :: 
the equation 

Substituting P - "P t -fi» and V* t-Jt^ into the first equation and 
subtracting out the steady state equation, we arrive at the result 



-p/CL - -pT = -IT* 

Next we consider the short section between (2) and (2 1 ). 
Neglecting the velocity of the fuel in the tank, we assume that 
Bernoulli's equation can be applied to the nonsteady as well as 
steady state flow when the perturbations are small, 

2- 

Pa = "Pa' 




BCL 1005 

Here 

9 Ug.' * ^ 2 y * flow velocity '"-'ft/sec. 

Subtracting and neglecting the second order of the perturbation -*-*-^» » 
we obtain 

To express this equation in terms of the weight flow rate W lb/.pec 
we substitute 

and arrive at the result 

-Pa * -PV + ( YV A? ) < 2 > 

Continuing from (2 1 ) to (1) and assuming a pressure loss due to 
friction in the form (L S£~ l7 ^ > Newton's second law yields 

•P, - "Pa/ " -fif * "^T TJ,** T^ft, 

where the acceleration of the fluid in the line is 
Bernoulli's equation for the steady state then becomes 
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Recognizing that ==■ O > and TXj = 9 the subtraction of 

the steady state equation yields 

or in terms of the weight flow rate 



- -fV T — ; " " ' ru ' <3) 



The quantities and , can now be elimiated from Equations (1) 

(2) and (3), and letting p. = p = suction pressure at the pump inlet, 

x s 

the resulting equation is 



if jar, f fa ♦ fy^y a = ^ T -p s (*) 

Its Laplace transform with zero initial conditions is w 



(5) 

It is worthy of mention at this point to call attention to the 
fact that the left side of Equation (5) is in terms of the inertiance 
and resistance defined as 

*— f\k A — inertiance 




R * 0 



resistance 
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The compliance in the above equation is zero since the fluid was assumed 
to be incompressible. 

In order to identify the transfer function for the block diagram, 
we rearrange Eq. (5) as follows 



La ( J 
which is represented by the block diagram of Figure 3. 



f-r 







V 
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>. 


i 
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Li ^S+-^ 





Figure 3- Block Diagram of Fuel Feed System 
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(b) Turbo Pump Including Cavitation Bubble 

The pump performance is limited by cavitation, or the formation of 
gaseous bubbles which occurs when the fluid pressure at any point becomes 
less than the vapor pressure of the fluid. Cavitation is most likely to 
occur at the pump inlet or the leading edge of impeller vanes. Bubbles 
collapse when the fluid pressure again reaches a value above the vapor 
pressure 4 Excessive cavitation will produce fluctuation in the flow rate 
and erratic combustion in the rocket engine which it supplies. 

Fuel entering the pump is accelerated towards the impeller periphery 
and enters the volute at high speed. The expansion in the volute then 
converts the large kinetic energy due to velocity into potential energy 
or pressure head c 

Pump characteristics are generally displayed in terms of pressure 
head rise against the total pressure head at the suction side and the 
flow rate. Typical curves are shown in Figure 4 and Figure 5. The total 



pressure head at the suction side is equal to (f 
and is called the net positive suction head NPSH> 




A 



X 




AH 



Figure 4 - Head versus Inlet Pressure 
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Figure 6 - Pump Flow Rate with Cavitation Bubble 




Figure 7 - Volume-Pressure Relation for Cavitation Bubble 
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Cavitation, or the formation of bubbles at the pump inlet,, increases 
as the NPSH decreases, Assuming the total volume of the bubbles to be "V^ 
we can apply the principle of continuity of weight flow rate between the 
$uction and discharge ends of the pump shown in Figure 6, 

Since the volume Y k is that of a large number of small bubbles which 
form and collapse in the inlet region, it is doubtful whether a polytropic 
gas law, such as * constant, can be applied.. However, the 

relationship between the volume *V U and the suction pressure 
must be of the general form shown in Figure 7* 

The cluster of bubbles in the region of the pump inlet can be viewed 
as a liquid-gas mixture of large compliance in comparison to that of any 
liquid itself. Thus, we have a column of liquid supported by this spongy 
mixture which behaves like a spring and it would be of interest to determine 
its spring rate ^ « 

The spring rate of a mechanical spring is determined by applying an 
increment in the load and measuring the deformation. Thus, if we assume 
the cross sectional area of this spongy mixture to be constant and equal 
to that of the pipe at that point, the spring rate can be. expressed as 

K b - —ET~ T% (9) 

If T§ is know, then the slope ..... * may be found from . . ? 



Figure 7, 



AY, 
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The spring rate \\ ^ may also be related to the natural frequency 
of the liquid column described above. Tho mase:;of the supported column 
is (JLji + and hence 



K 



A value of K corresponding to the first natural frequency of the 
missile would certainly represent an unfavorable compliance to the 



system* 



To complete the transfer function for the pump including the effect 
of the cavitation bubbles we obtain from Eq, (9) 

Eq* (8) can then be written as 

^r d * ^ - do) 

Taking the Laplace transform of Eqs. (7) and (10) the transfer functions 
for the pump and the cavitation become 




(id 



^4 



^c-b = > 5 ^ - -j^-s^ct^ (12) 



The block diagram representing these equations is shown in Figure 8. 
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Figure 8 - Block Diagram of Turbo-Pump Including Cavitation Bubble 



(c) Discharge Line and Injector 



The discharge line carries the fluid under high pressure from the 
turbo-pump to the injector and combustion chamber, as shown in Figure 9 




Figure 9 - Discharge Line and Injector 
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Applying Newton's second law to the -discharge line from d to d 1 , 
assuming the pressure loss due to friction to be C i -yrr L) , A 

We subtract from the above equation the steady state equation 

o = P d - % - Ca^u/ 

and introducing the relationship 



we obtain the result 



— Jfr = ^"-f^, - (13) 



Proceeding to the injector we write Bernoulli's equation between 
d' and the combustion chamber where the pressure is T~£, 

■ 
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Again subtracting the steady state equation and replacing (J in terms 
of "NaT* we arrive at the equation 



(14) 



Eliminating ' between Eqs , (13) and (14) and taking the Laplace 

transform, the transfer function for the discharge line and injector 
becomes 



'A 



A2£ - A 



(15) 



Its block diagram can then be drawn as in Figure 10. 



-ft. 




Figure 10 - Block Diagram of Discharge Line and Injector 
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The ignited gas in the chamber travels through the converging- 
diverging nozzle which is designed so that sonic speed is reached at 
its throat. The expansion then takes place supersonically to the exit 
with the developed thrust F~ * This process from chamber pressure ~Pq, 
to thrust F" is again dependent on many factors and hence it is usual 
practice to simplify the result in terms of a thrust coefficient C f 
which is a function of the specific heat ratio, the ratio of exit to 
throat area:;f\p y^^and the ratio of exit to chamber pressure T^jT^ . 
The equation is then written as 



J" ~ 



and its transfer function is 



- = C F A t 



(18) 



The block diagram for the engine can then be draw as in Figure 11, 





C 


JVpS 









— • — 




Figure 11 - Block Diagram of Rocket Engine 
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(e) Structural Dynamics 

The dynamical equations for the structure can be established from 
Figure 2, Since we are interested only in the fluctuation due to the 
variation in thrust, these equations become 



Taking the Laplace transform with zero initial conditions 



Thus the transfer functions for the structure become 



(19) 



where 



o 



Its block diagram can then appear as in Figure 12, 
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Figure 12 « Block Diagram of Structure 



In a more sophisticated approach, the above simplified technique 
can be replaced by the mode- summation method utilizing generalized 
coordinates. The deflection at any point is then represented as 



where 



and 



= normal mode of the structure 



Structural dynamics would still be represented by a block diagram where 
the thrust p is the input and the displacements at 1 and 2 are the output, 
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System Block Diagram 

It is a simple matter now to assemble the various block diagrams 
into a closed loop circuit shown in Figure 13. System behavior can then 
be computed by programming on the analog or digital computer. 

Computation and Evaluation 

Numerical values for the system parameters must be established for 
the duration of flight of interest to this study. Changes :'.in:these 
parameter values iaay be incorporated in the computation to establish 
their effect on the behavior of the system, 

"Fix-devices" which have been used successfully are the piston- spring 
accumulator and the standpipe accumulator at the suction end of the fuel 
feed line shown in Figure 1. Such devices may be incorporated into the 
system diagram to study their effect on performance. The increase in the 
gas pressure ~P-y in the tank will also help to stabilize the system, in 
that it will increase the pressure at the pump inlet. and reduce cavitation. 



Approximate Numerical Values 

I 40 ~ - , 

L j - 0.3 
j^L & r 1900 

^ A.p ' 90 x 32.2 x 0.40 2 



^ 2.2 



*-^[^\ 0,2 
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- 0 



0L+ l ^ 2 

13 p = 90 

CV 0 V 134 

* 0,0015 ' 

^ * 0 

= "rf* ^ ~ ■ 76 for unfavorable condition of 12 cps. 

^b = 762 x 32?2 x 44 * 710 > 000 lb / ffc2 
-2, 



0.0000204 



90 x ,40 2 <n> 

K b " ?io,ooo 

Problems for Students 

(1) With the bulk modulus of the fluid defined as B = , show that 

the compliance of the fluid in the pipe of area A is * What is 

B 

the compliance per unit area for water? 

(2) Derive the transfer function for the standpipe accumulator and indicate 
how it is fitted into the block diagram of Figure 13 ♦ 

(3) Using the mode summation technique, determine the transfer functions 
for the structural dynamics and draw its block diagram to replace the 
one in Figure 13, 
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Observed Natural Frequenc 

Nat, Frdq, 105 
rad/sec 65 

Cd 2 rad/sec 119 
-V^ c 60,000# 
V 2 " 90,000# @ t 
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Time-Sec 

115 125 135 145 



68 74 83 105 - 7 (- 

136 159 173 174 

constant 
■ 105 sec. 
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Notes for Instructor 

The mode summation method assumes the displacement at any point to 
be expressible as 

where is the i*"^ normal mode shape and ^i^"^-) ^ s the generalized 

coordinate for the i*^ mode satisfying the equation 

In the^above e quation is the natural frequency, ^> the damping factor, 

the generalized mass, all for the i** 1 mode, and |- (-^ 



c 

is the applied force ♦ 



For the Pogo problem is applied at the end of the missile 

~ O , and the generalized force becomes 

A 

o 

The Laplace transform of the displacement is 

where 

Thus at 1^*^. and we have ^ n A > A 

<£ filling! £^±L> 
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Since data are given for only the first two modes the result with zero 
damping becomes 



The block diagram then becomes 























(P-2. <o> UP?- 


















(P, Co^> 































^1 



which replaces Figure 12. 
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Equations for the Standpipe 



P A V 



> >»<i + 



Newton's second law applied to the fluid in the standpipe is 
Eliminating the steady- state equation we are left with 



X L Z - ^ - -fV - ^ ^ 

The trapped air in the standpipe must obey the gas law 



(i) 



Expanding by the bionomial theorem and neglecting second order 
terms, we have (V^f-J^^^ ~K 7^%^ + - ^ = Po-V^ 



(2) 
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Eq. (1) now becomes 



We next introduce continuity of flow 



The transfer functions from Eq, (3) and (4) are 



(3) 



(4) 
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It should now be recognized that we can let a* "P^k and that 
JjJ^ in Figure 3 must be changed to -U/| . Thus the block diagram for 
the standpipe becomes 



A" 



fAS 


1 




1 







^5 



K 



b i 



.— i 



and its position in the large block diagram is also indicated. 
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